arXiv: 1502.03076vl [hep-th] 10 Feb 2015 


On thermal fluctuations and the generating functional in 
relativistic hydrodynamics 

Michael Harder, Pavel Kovtun, and Adam Ritz 
Department of Physics and Astronomy, University of Victoria, Victoria, BC, V8P 5C2, Canada 

(February 2015) 


Abstract 

We discuss a real-time generating functional for correlation functions in dissipative rela¬ 
tivistic hydrodynamics which takes into account thermal fluctuations of the hydrodynamic 
variables. Starting from the known form of these correlation functions in the linearized 
regime, we integrate to find a generating functional which we can interpret within the 
CTP formalism, provided the space-time and internal global symmetries are realized in a 
specific manner in the (r, a) sectors. We then verify that this symmetry realization, when 
implemented in an effective action for hydrodynamic fields in the (r, a) basis, leads to a 
consistent derivative expansion for the constitutive relations at the nonlinear level, modulo 
constraints associated with the existence of an equilibrium state. 

1 Introduction 

Hydrodynamics is a low-energy effective description of many physical systems in local ther¬ 
mal equilibrium. Relativistic hydrodynamics for normal fluids is a set of partial differential 
equations expressing the conservation of the energy-momentum tensor and other conserved 
currents, such as the baryon number current. Such a description is classical, rather than field- 
theoretic, and one can enquire about the low-energy effective field theory corresponding to the 
hydrodynamic regime, and how it is related to the classical hydrodynamic equations. A funda¬ 
mental object in field theory is the generating functional W[A, g] where and = Vnu+hfiu 
are the external sources (gauge field and the metric). The variations of hF[A, 5 r] with respect 
to the sources give rise to hydrodynamic correlation functions, i.e. to correlation functions 
of the energy-momentum tensor and the current in the limit of small frequency and 
momentum. 

The motivation for finding such a VF[A, 5 r] is the following. There is a plethora of n-point 
real-time response functions, differing by time ordering and symmetrization of the correspond¬ 
ing operators. These correlation functions can be conveniently classified in the closed time 
path (CTP) formalism by labeling the operators according to the two parts of the time con¬ 
tour [1, 2]. The response functions computed in classical hydrodynamics by varying the 
on-shell and with respect to the sources are the fully retarded functions, or raa...a 
functions, in the notation of Ref. [2]. While for two-point functions, all response functions can 
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be reconstructed from the ra function by the fluctuation-dissipation theorem, the same is not 
true for general n-point functions [2], In classical hydrodynamics, the fluctuation-dissipation 
relations are not contained in the hydrodynamic equations, and have to be imposed by hand. 

Second, classical hydrodynamics misses information about thermal fluctuations of the hy¬ 
drodynamic degrees of freedom themselves. Such real-time thermal fluctuations will lead to 
effects (long-time tails and momentum-space non-analyticities in two-point functions) which 
can not be captured by classical hydrodynamic equations [3, 4]. Relativistic fluids are not 
immune to such effects, though the fluctuation corrections become suppressed in the large- 
limit [5, 6]. A full real-time generating functional must capture both the fluctuation-dissipation 
theorem, and the effects of thermal fluctuations. 

One way to arrive at a hydrodynamic effective action is to modify the hydrodynamic 
constitutive relations by introducing random stresses and currents whose correlation properties 
are chosen so that the fluctuation-dissipation theorem is satisfied in equilibrium [7]. Such a 
construction is phenomenological: it takes the hydrodynamic equations as given, while the 
random stresses and currents become extra dynamical degrees of freedom, to be integrated 
over in the hydrodynamic path integral. It is not immediately clear how to proceed with the 
systematic derivative expansion and the coupling to external sources in this formalism. See 
e.g. [8, 9] for recent discussions in the context of relativistic hydrodynamics. From a field- 
theoretic perspective, it would be more natural to implement the hydrodynamic equations, the 
derivative expansion and the coupling to external sources at the level of the effective action 
which respects the relevant symmetries of the microscopic physical system. 

The present paper is a step in this direction. We start with the bottom-up approach 
in Section 2, asking a simple question: what is the generating functional that gives rise to 
the known hydrodynamic two-point functions of linearized hydrodynamics in equilibrium? 
Incorporating the appropriate background sources, a structure emerges that is consistent with 
expectations from the CTP formalism, expressed in the so-called (r, a) basis. Going beyond 
linearized hydrodynamics requires understanding the symmetries of the effective action. In 
the CTP formalism, the set of symmetries is doubled (call these symmetries and Q 2 ), 
corresponding to the two branches of the time contour. The classical hydrodynamic equations 
on the other hand manifest only one symmetry, the diagonal (physical) Qr- We discuss the 
symmetries and derive the relevant Ward identities in Section 3, which generalize the classical 
hydrodynamic conservation equations. These classical equations have the schematic form of 
conservation laws 

V-Jr = 0, (1.1) 

where Sr = , etc. and we have ignored possible explicit symmetry breaking terms that 

may be present on the right hand side. Based on the results of the bottom-up analysis, we 
further argue that the hydrodynamic effective action can be built from the degrees of freedom 
that arise in a low energy nonlinear symmetry realization, analogous to the spontaneous 
breaking of Qi x Q 2 ^ Gr- The degrees of freedom of the effective theory thus include modes 
analogous to the extra Goldstone modes arising from symmetry breaking. The effective action 
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has the following schematic expansion in terms of the a-sector (fluctuation) fields (/Jq, 

'S’efF =^r •Jr ' + . . . . (1-2) 

Here Jr and K-r are functionals of the r-sector (physical) fields, with indices suppressed. 
'Dipa denotes the appropriate covariant derivative of <pa, which is manifestly invariant under 
the hidden symmetries Qa, orthogonal to Qr- The term linear in Dipa, on varying with respect 
to ipa, enforces the classical equations of motion D ■ Jr = 0, a generic feature of the CTP 
action. The term quadratic in Dipa provides the additional structure necessary to satisfy the 
fluctuation-dissipation theorem, and can be interpreted in terms of supplying fluctuations. It 
is apparent that the derivative coupling of ipa is consistent with the dependence expected for 
Goldstone modes associated with a nonlinear symmetry realization. We test this expectation 
at the nonlinear level by writing the general CTP effective action in a low energy derivative 
expansion for both the classical hydrodynamic fields, and for the fluctuation modes ipa- In 
Section 4 we implement the expansion to first order and identify all the expected transport 
coefficients of first-order hydrodynamics in terms of the parameters of the effective action, 
modulo certain constraints associated with the existence of an equilibrium state. We conclude 
in Section 5 with a list of open questions that need to be resolved in order to have a complete 
picture of the hydrodynamic generating functional. 


2 Bottom-up approach 
2.1 Diffusive mode 


We start with linearized hydrodynamics, and consider the simplest hydrodynamic process 
which is diffusion. It is described by the diffusion equation 


dtn — DV^n = 0 , 


( 2 . 1 ) 


where n is the charge density fluctuation, D is the diffusion constant, and = did^. The 
linear response theory gives the following two-point functions of the charge density n{t, x) in 
thermal equilibrium: 


Gra (^) h) 


iuj — ’ 


Gariys k) 


iuj + D]J ’ 


Grr h) 


—4i TDx 
a;2 (Dk2)2 ' 


( 2 . 2 ) 


The first one is retarded, the second one is advanced, the third one is {—i) times the anti¬ 
commutator, and Gaa is identically zero. Here T is the equilibrium temperature, and x = 
{dn/dfi)fj,=o is the static charge susceptibility. What is the effective action which gives rise to 
these correlation functions? 

In relativistic hydrodynamics, the diffusion equation emerges from the current conservation 
equation = 0 in the Landau-Lifshitz frame [10], linearized in small fluctuations close to 
the equilibrium state at zero chemical potential. We can couple the system to an external 
gauge field which gives rise to </^[^]. The variation of the hydrodynamic on-shell current 
with respect to the source A can give rise to Gra and Gar, but not to Grr- This is because 
the linearized current conservation equation in the presence of the source has in it D and 
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a = Dx, but not T. In order to get Grr, hydrodynamic equations coupled to sources are not 
enough, and one has to use extra information, namely the fluctuation-dissipation theorem. 
The effective action must incorporate both the hydrodynamic equations with sources, and the 
fluctuation-dissipation theorem. 

It is intuitively clear that there is no local path integral action S[(j)r,(l)a\ quadratic in the 
fields, which would give Gra = —i{4‘r4‘a), Gar = —i{4‘a4‘r), Grr = —2i{(j)r4>r) with the response 
functions (2.2). This is because the response functions (2.2) have in the numerator, hence 
the corresponding action will have in the denominator, which is not local in space. For a 
quadratic action of the form 



in the path integral, the matrix 

- T)k2 2iT ) 

gives the correct response functions (2.2). Here, the upper left element is rr, upper right is 
ra, bottom left is ar, and bottom right is aa. In particular, 

{(j)a4^0) i{P ')a0 ) 


where the indices run over r, a. The action (2.3), (2.4) is not real: it is complex, but in such 
a way that the functional integral with the weight converges. Clearly, the matrix (2.4) is 
not analytic in k, and the quadratic action S[(j)r,<pa] is not local in space. 

The action can be made local by introducing auxiliary fields. Let us define a new field (fa 
as (j)a = Dx^‘^^a- Consider the following action 


5[0,.,0a,¥5a, A] = dtd'^X (faidt - D'V‘^)(f)r - iipaT(j)a +Xi(j)a - Dx'V'^^a) 


(2.5) 


where the auxiliary field A is used to impose the constraint which defines ifa- We can further 
define the generating functional as 

Z[ar,aa] = , 


where the effective action is given by Eq. (2.5), and Ur, Oa are external sources. The physical 
meaning of cj)r is the density fluctuation field n. By construction, this generating functional 
reproduces (2.2). Integrating out A and (pa leaves 




I DprDifa exp 


i / [(fa {dt(pr - D'V'^cpr Dx'V'^ar) 
• Jt,x. 

-i^a TDxV^ <Fa + CLaP^r] ■ 


( 2 . 6 ) 


The source now appears as a correction to the equation of motion of cpr- This is precisely 
how the source term should appear, based on the full hydrodynamic description. Indeed, 
in relativistic hydrodynamics, the diffusion equation (2.1) arises from current conservation. 
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To first order in the derivative expansion, the constitutive relation for the current in the 
Landau-Lifshitz frame is 


, (2.7) 

where is the fluid velocity satisfying = —1, T is the local temperature, /r is the 

chemical potential, and + u^vy. The electric held is , where F^y is the 

held strength. Let us turn on only. For linearized huctuations around the equilibrium state 
with = (1, 0), T = const, and /r = 0, we have = n, = —ad'^fj. + ad'^Ao. Charge density 
and chemical potential huctuations are related by n = y/r, and the electrical conductivity is 
a = Dx- Current conservation = 0 now gives 

dtn — D'V^n + = 0 , 

modifying the diffusion equation by the source term proportional to V^Aq. Thus, by com¬ 
paring with (2.6), we identify with Aq. The hydrodynamic equation coupled to the source 
contains both D and y, but not T. The dependence on temperature comes from the kinetic 
term for ipa in the effective action (2.6). 

One can rewrite the kinetic term for ipa as 

This makes the action linear in (/?„, enforcing the equation of motion for (fir with the Gaussian 
noise in the right-hand side. Integrating out (pa gives 

where n[ar,r] in the exponent is a solution to dtn — D'V^n + Dx^'^cir = —d^rk, for a given 
source ar{t,x) and noise prohle rj(t,x), and the average is over the Gaussian noise r^. This is 
the standard relation between stochastic equations and path integrals. 

One can write down a covariant generalization of the generating functional (2.6), 

Z[A^, Aa] = jo^rDipa ArA,] ^ (2.9) 

where 

S = Jdtd‘^x{J^^[4>r,Ar]D^ipa + iTaAf^''Df^(paDy(pa) , (2A0) 

and Dfj^(pa = d^ipa + A“. The gauge held A^ is the physical gauge held, while A^ is only 
used as a tool to access correlation functions, and needs to be set to zero at the end of the 
calculation. The effective action is invariant with respect to both r-type and o-type gauge 
transformations. The conserved current obtained by varying the effective action with respect 
to Aa is the classical hydrodynamic current plus the huctuation correction, 

+ 2iTaA>^'^Dyipa . ( 2 . 11 ) 


= / Dr e^TDx 


friri^ijipadkrk 


( 2 . 8 ) 
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2.2 Shear and sound modes 


Let us now look at the hydrodynamic shear and sound modes. It will be easiest to work with 
an uncharged fluid. In the Landau-Lifshitz frame, the energy density fluctuation e = — e 

(with e the equilibrium value) and the momentum density fluctuation vr® = obey 


dte + dfcVTfc = 0 , 

dm + Vgdie - MijTTj = 0 , 


where = dpjde is the speed of sound squared, Mjj = 7 ^(V^( 5 jj — didj) + 'jgdidj, and the 
damping coefficients are 7 ^ = r]/{e + p), 7 ^ = (/{e + p), 7 <i = 7 (^ + The hydrodynamic 

retarded functions are 


^ ra 

^TTiTTj 


(w,k) = ( 


WPv 


k2 


tOJ 


-7r?k2 


+ e 


+ 


kikj 


WUJ^ 


1^2 _ \p2y2 _j_ 


P 


( 2 . 12 a) 


G'™(u;,k) 

G™(a;,k) 


G:“(tu,k) = 


wuki 


^2 _ ]^2^2 _j_ ’ 




^2 _ ]^2y2 _|_ ■lijj^^\p2, 


- e, 


( 2 . 12 b) 

( 2 . 12 c) 


where w = e + p. These functions are obtained by first solving the hydrodynamic equations in 
the external metric, and then by varying the resulting solution for with respect to 

the metric, see for example [11]. Note that (u;,k) is analytic as k —)• 0. All aa functions 
vanish identically, while the equilibrium rr functions can be obtained from the fluctuation- 
dissipation theorem as 

4?T 

G7b = —ImG7s. (2.13) 

(jJ 

We would like to find an effective action which reproduces the above response functions, as well 
as the corresponding rr functions. Let us choose e and vrj as our variables. By analogy with 
the diffusive generating functional ( 2 . 6 ), one can make a guess for the generating functional 
for the shear and sound modes: 


Z[hr, ha] = [DeDTTiDi^QDipi exp i f [ipf [dm + v^die - MijiTj + w dth^i - dih-Qo) 

J ^ J t,X 

+ Po{dt(^+dk'K^) — iTw + |/ioo(e + e + ^e^oo) + • (2-14) 


The fields e, tTj here are r-type fields, and the auxiliary helds (/ 9 q, c/jf are a-type helds. The only 
sources turned on are Hq^. The /ig^ sources in the equation of motion come from = 0 

in the Landau-Lifshitz frame. The /ig^ sources come from Tr'^h^^. 

An exercise with Gaussian integrals, given in Appendix A, shows that the generating 
functional (2.14) does indeed reproduce the equilibrium response functions (2.12) of linearized 
hydrodynamics in an uncharged relativistic fluid. The fluctuation-dissipation theorem (2.13) 
and the vanishing of all aa functions automatically follow from the structure of the effective 
action in (2.14). Just as in the example of diffusion, the generating functional can be cast into 
the form of a stochastic equation with Gaussian noise. 
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One can write down a covariant generalization of the generating functional (2.14), 

Z[hr,ha\ = , (2.15) 

where 

S = . (2.16) 

Here ~ ~ To first order in the derivative expansion, the classical 

energy-momentum tensor is given by the standard expression in the Landau-Lifshitz frame, 

= eu^u'^ + , 

with + A^"A^^ - + (A'^^'A^^, and + u'^u’^. The 

indices are raised using gr, which is the physical metric, while ha is only used as a tool to 
access correlation functions, and needs to be set to zero at the end of the calculation. The 
effective action is invariant with respect to both r-type and a-type diffeo transformations. 
The latter act as position-dependent shifts of (p^. The conserved energy-momentum tensor 
obtained by varying the effective action with respect to is the classical hydrodynamic 
plus the fluctuation correction, 

^ 2iTG^^^P^Vpipl . (2.17) 

One can explicitly check that the generating functional (2.15) gives the correct equilibrium two- 
point correlation functions for all components of the energy-momentum tensor in linearized 
relativistic hydrodynamics, and that the fluctuation-dissipation relations are satished. This 
is a non-trivial check of the validity of the effective action (2.16) for linearized hydrodynamic 
fluctuations. 

The structure of the effective action (2.16) is easy to discern. Suppose the kinetic terms 
for were not there. Then integrating over the would impose = 0 as an exact 

operator equation. The generating functional then becomes 

Z[hr,ha] = 

where stands for T'^'^[T[gr],u[gr], gr]- Taking the variation with respect to 

produces one r-insertion of the on-shell , and subsequent variations with respect to h'" 
will produce raa ... a hydrodynamic functions, obtained in the standard way by varying the 
on-shell energy-momentum tensor. The kinetic terms for 99 “ are responsible for the fluctuation- 
dissipation theorem, and allow one to evaluate correlation functions with more than one r- 
insertion. These terms are responsible for thermal fluctuations of the hydrodynamic modes 
allowing them to go off-shell, and will give rise to hydrodynamic loop corrections and running 
of transport coefficients. 

3 Top-down approach 

In order to correctly extend the generating functional W = —ilnZ beyond Eq. (2.15) of 
linear hydrodynamics, one needs to be more systematic about the underlying symmetries. 
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The effective action identified in Section 2 has the schematic form, 

'S'eff =^r •Jr ' + • • • (3-1) 

Here Jr and fCr are functionals of the r-sector fields, with indices suppressed. The term X^, 
which is independent of (fa, was not present in the earlier linearized discussion, but we include 
it here for completeness of the expansion and will return to it below in the context of the path 
integral measure. The term linear in T>ipa, on varying with respect to ^pa, enforces the classical 
equations of motion V ■ Jr = The term quadratic in Vipa supplies fluctuation corrections, 
and provides the additional structure necessary to satisfy the fluctuation-dissipation theorem. 

3.1 Non-equilibrium CTP contour 

We identify the expansion (3.1) as one that appears naturally within the nonequilibrium 
Schwinger-Keldysh CTP formalism, involving a doubled set of fields and symmetries corre¬ 
sponding to two time contours, see e.g. [1, 2]. 

For a quantum-mechanical system with fundamental degrees of freedom q, which at time 
to is characterized by the density operator p, the CTP generating functional is given by the 
path integral over the fundamental fields 

Zp[ji-,j2] = j dqidq2dqf {qi\p\q2) j DqiDq2 To L{q2,j2) ^ 

where ji and j 2 are the external non-dynamical sources, and the boundary conditions are 
= qi, <72(^0) = Q 2 , = Q 2 {tf) = qf- The generating functional satisfies Z\ji,ji] = 1 

as well as Z[ji,j2]* = Z[j2,ji], thanks to trp = 1 and p = pT In terms of the (r, o) variables 
Qr = {qi+q 2 )/‘ 2 , qa = qi-q 2 , the action is 

-S[q2,j2] = jqaE{qrJr) + 0{ja,ql), (3.3) 

where E{qr,jr) is the classical equation of motion, see e.g. [12, 13]. For a thermal equilibrium 
state, the temperature dependence comes from the matrix element of the density operator. A 
symmetry Q of the classical action will lead to a doubled symmetry Qi x Q 2 of the generating 
functional provided the integration measure is invariant, and the density operator p transforms 
covariantly. For the latter to be true, Gi x Q 2 must reduce to the diagonal Gr at the initial 
time to) as there is only one (physical) symmetry characterizing the initial state.^ For local 
symmetries characterized by a continuous parameter we have 

J2] = Zp>[ji -F J2 + d^ 2 h] , 

with ^i(to) = 6(to), ?i(t/) = 6(t/), and p' = 5^po)p. 

We seek to find a similar two-source generating functional for low-energy excitations of 
near-equilibrium states, in which case the effective degrees of freedom are the hydrodynamic 

^Note that the doubled symmetry of the generating functional does not mean that the symmetries of the 
theory magically double. The two-source generating functional is just a means for convenient classification of 
correlation functions, and the time contour can be chosen to run back and forth more than once. 
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Figure 1: A contour with two times and two /3’s. 

modes (in place of qr) and the corresponding auxiliary helds (in place of qa)- The classical 
equations of motion for the r-type fields are the conservation laws for the energy-momentum 
tensor and other currents. Note that the condition Z[ji,ji] = Ij ensured in the microscopic 
theory by the normalization of the initial density matrix, may result in a nontrivial measure 
factor at 0{a^) when we consider a low energy effective action in place of the microscopic 
action (3.3). In this way, we observe that the a-sector expansion of (3.1) can naturally emerge 
from the CTP formalism. The doubled symmetry Qi x Q 2 of the microscopic description (with 
the off-diagonal Qa broken by the initial state) needs to be realized in terms of the effective 
degrees of freedom. We will subsequently propose that the correct symmetry realization in 
the hydrodynamic effective theory is a nonlinear realization of Q\ x Q 2 with an explicit (linear) 
realization of the diagonal Qr- The auxiliary a-type degrees of freedom then couple in the 
manner expected for the corresponding Goldstone modes. 

3.2 CTP metric sources and conservation laws 

To analyze the hydrodynamic regime, we need to consider the CTP formalism in the presence 
of sources for the charge current and energy momentum tensor. To that end we hrst introduce 
two metric sources and so that the generating functional is 1F[5^,5^]. We demand 
that W is invariant under two sets of diffeomorphisms: Pi which only transforms 5^, and T >2 
which only transforms . We anticipate that there will be additional fields in the theory, such 
that at low energies we can identify physical (r-sector) fluctuating modes associated with the 
hydrodynamic degrees of freedom, e.g. the temperature, fluid velocity, etc. We will not need 
to specify these modes explicitly, but we necessarily assume that their dynamics is consistent 
with T>i and T> 2 . As one example, we might envisage a system with an initial equilibrium 
state, with inverse temperature j3i and ^2 in the two sectors. The corresponding time contour 
is depicted in Fig. 1. 

We will be working with linear combinations of the sources which provide easy access to 
the retarded and symmetrized functions. We define 

9 — 2 [9 + 9 ) j 9 — 9 — 9 ■ 

The g'' source corresponds to the physical metric. The ^7“ source needs to be set to zero 
at the end of the calculation. Diffeo invariance of W gives rise to conservations laws of the 
energy-momentum tensor. We define 

5gW\g\g^] = I ^^y^{TnSgU " / , ( 3 - 4 ) 
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where 


= gU^uCi + + ^^gU , ^glu = glxdui2 + gl\^^^^2 + dxg% C2 

are the variations of the metric. Diffeo invariance of W gives 

v^'(rn = o, v2(rr) = o. ( 3 . 5 ) 

These conservation equations can be expressed in the r and a basis. To that end, we define 


{TD = i {TD + 5 {T2n , {TD = (^D - (^ 2 '^") • 

The average {Tr'^) is the physical stress tensor. Thus 

5gW[g\g^] = I{TDK-^ + J■ (3.6) 


The variations Sg'^ and Sg^" can be expressed in terms of = (?i+? 2)/2 and = Ci—■^ 2 . For 
the r-metric we have 5g'^^^ = drg'^^y + ^agliy, where 

^rg^iu = g^iX^i'^r + glx^g-^r + dxg^iu Cr ) ^ {g^gX^u^a + g^X^fiCa + dxg'j^u ■Ca) • 

Similarly, for the a-type metric we have <5^“ = drg^u + ^aguu where 


^rg‘liu — g'^X^i'ir + g^X^fiCr + dxg'^y C. 


a 

r 1 


^ag% — g^iX^y^a + glx^g^a + dxgliu Ca • 


Diffeo invariance of the generating functional (3.6) then gives rise to the conservation laws for 
{T^'^) and (T^^). Upon setting 5 “ = 0 these reduce to 


v^(rr) = 0 , v^(rr) = 0 , (3.7) 

where by we denote the physical covariant derivative, evaluated with respect to g^. 

It is straightforward to incorporate external gauge fields and the corresponding 

combinations = \{A + A), A°' = A — . The currents are dehned by 

5aW = j^{j^)5Al- jx/^{4)5Al. 

Gauge invariance in the 1 and 2 sectors leads to 

v;,(jr) = o, vI{4) = q, 

while diffeo invariance leads to (3.5) with the usual Joule heating terms in the right-hand side. 
Again, the conservation laws can be expressed in the (r, o) basis. We define 

{4) ^ {4) = v^(Jf) - xf^44) • 

The average {Jr) is the physical current. Gauge invariance of IT then gives 


v^(jn = o, v^(jn = o. 
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The diffeo transformation properties of the r and a gauge fields are 6 A''^ = with 


6rA^=Cd,A^.+Ald^C, 


SaA^^ = l{CduA';, + A:d^C) , 


as well as 5^“ = 6 rA^ + SaA^j^, 


with 


6rAl = Cd,Al + Ald^C , SaAl = CduAl + Ald.Ca 


(3.8a) 


(3.8b) 


Diffeo invariance of the generating functional then gives rise to the following conservation laws: 


g:,v,{Tn + i {gtxd.iTD++ 5 “Ar^;rr) = F^AJn + kp^AJa ), (3.9a) 

alxSATn + aAdATD + n.uiTD + 5“Ar^;(rr) = FxA^n + > (3-9b) 

where Tx/iu = ^(^/.jfl'i/A + dyg^iX — dxgi^u)- We have written the conservation laws in this form 
to avoid using the inverse of the a-type metric. Taking further variations with respect to the 
metric will give rise to Ward identities for two- and higher-point correlation functions in the 
ra basis. 


4 Derivative expansion for the effective action 

To map the general analysis of the CTP effective action onto classical hydrodynamics, we need 
to understand how diffeomorphism (and/or gauge) invariance is realized in the low energy 
(hydrodynamic) regime. 

Let us first consider extending the straightforward realization of diffeomorphism invariance 
in an equilibrium state [14, 15] to the hydrodynamic regime. We can characterize all near 
equilibrium states in terms of a timeline vector field /3^ = where ,5^ is a timelike 

Killing vector characterizing the equilibrium state. If we denote the generating functional in 
this state as r[( 7 , /3],^ then for sufficiently well-behaved states, a derivative expansion for r[gf, /3] 
exhibiting manifest diffeomorphism invariance can be implemented. As an example, consider 
terms up to first order in the derivative expansion. We take /3^ = u^/T, with T = l/\/— /3-/3, 
then to first order 

T = j(A+^x {p{T) + a(r)V-u + b{T)t + ...), 

where T = u^dfj^T, and p, a, b are arbitrary functions of T. In fact, integration by parts shows 
that the two structures T and V • u are not independent; the independent coefficient can be 
identified as the combination {b — a'). Varying T with respect to the metric (keeping /3^ fixed) 
gives the following energy-momentum tensor: 

Tt^u ^ ^ + (5_a') , (4.1) 

frame—dependent 

^In equilibrium states, the generating functional of [15] is identified here as r[y, P] = W[g, S'] — fy/—g /3^S^, 
where S'' is a source for the held so that Por equilibrium states with no source, S)_i = 0, 

the energy-momentum tensor is conserved. 
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where + u^u'^ is the usual transverse projector. Identifying p{T) with pressure, 

the hrst two terms give the standard ideal hydrodynamics, as one can see from the thermody¬ 
namic relation Tp' = Ts = e + p. The term proportional to [b—a') looks like a bulk viscosity 
contribution, however it is an artefact of using the thermodynamic frame: the frame-invariant 
one-derivative scalar [16] vanishes, and the {b—a') term gives no bulk viscosity and no dissi¬ 
pation. Indeed this symmetry realization appears essentially Euclidean, and in equilibrium 
states gives Euclidean (zero frequency) correlators which analytically continue to retarded 
correlators in real time. 

In order to describe dissipative hydrodynamics, we need to add the a-type sources in the 
right way and some dynamical fields to the theory, in order to produce correlation functions 
which are not polynomial in spatial momenta. As noted earlier, the bottom-up construction 
(2.10) and (2.16) suggests a specific realization of the doubled symmetry Vi x P 2 in the 
hydrodynamic effective theory. Specifically, we assume an a nonlinear realization with an 
explicit (linear) realization of the diagonal T>r analogous to the treatment above. The presence 
of non-dynamical metric sources implies that the global symmetries of the theory are in effect 
weakly gauged. 

We can be agnostic about the precise mechanism via which this symmetry realization 
arises from the microscopic theory, and carry out a general parametrization of the low energy 
degrees of freedom. This is the approach we take in this section, making use of the natural 
derivative expansion in the hydrodynamic regime. See [17, 18] for a discussion of coupling the 
hydrodynamic degrees of freedom to a different type of Goldstone modes. 

4.1 Charge current 

Before considering diffeomorphisms in detail, to gain some intuition for the framework outlined 
above we first consider the simpler case of U{1) charge diffusion in flat space. The above 
picture suggests that we should consider t/(l)i x U{ 1)2 —)■ U{l)r in the hydrodynamic regime. 
Denoting the corresponding ‘Goldstone-like’ mode on general grounds the low energy 
effective action will have the form 


5eff[A^, Pa] = 5eff[e^, (4.2) 

where = D^pa = d^Pa + Aq is gauge invariant under the off-diagonal U{l)a symmetry, 
while Fr'^ = d^A’^ — d’^Ar is gauge invariant under the residual U{l)r- See Ref. [19] for a 
(technically) similar approach to the effective action for zero-temperature superfluids. The 
derivative expansion for the Goldstone modes is then equivalent to the expansion in the a- 
type fields. Note that consideration of any microscopic example of U{l)i x U{1)2 symmetry 
breaking provides some justification for the expectation that while ^eff [A,., ^a] rnay naturally 
be represented in terms of hydrodynamic degrees of freedom at low energy, this need not be 
the case if the effective action is written in terms of the 1- and 2-type degrees of freedom. 

To linear order in the a-fields we have 

5eff = y’j,^(A“+a^(^“) + 0(a2), (4.3) 
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where Jr = 63^^/5A^. The equation of motion for is then (9^J^ = 0 as required. In 
writing (4.3), we have dropped the O(a^) terms which only depend on the r-type fields and 
the r-type sources of the effective theory. These terms must ensure that the normalization 
of the generating functional ZlA^jA^'] = 1 for A^ = d^X with A(to) = 0 is preserved in the 
effective description. 

To determine the hydrodynamic constitutive relation for Jr, we write down a generic 
effective action depending on the allowed variables in the r and a sectors. We choose the 
variables for constructing the effective action to be given by the set {T,u^, n,^a, Fr'^}, where 
for the moment we assume a flat background geometry, T denotes a scalar that we identify 
with temperature, and is a normalized fluid velocity vector satisfying = — 1. We can 
identify the timelike vector (3^ = + (3'^ discussed above as f3^ = u^/T. In equilibrium, it is 

natural to instead start with (3^, which is then necessarily proportional to a timelike Killing 
vector, and construct the fluid velocity = (3 ^/^ and temperature T = 1/ ^/—/3^. This 
implies a particular dependence on the background metric as a source, which we will not need 
below. We choose to start with the conventional hydrodynamic variables T and for the 
following out-of-equilibrium analysis. 

In terms of this data, up to 0{a) there are three gauge invariant scalars at zeroth order in 
derivatives: 

ai=T, a2 = Ai, as = /^a , (4.4) 

where /Za = The identification of the a-sector chemical potential /Za with the invariant 

(sometimes referred to as the Josephson relation) reflects a redundancy in the set of 
invariants, as both determine the source dual to the conserved charge, see e.g. [20]. We will 
adopt the identification above to remove this redundancy. Note that ai,a 2 rsj 0(a°), and 
as ~ 0{a). Thus we find at this order, 

5® = jF{T,t,r,fia) + 0{a^) = I F^,u^I4^(^“ + 0(a2), (4.5) 

where F is the effective Lagrangian, and the derivative = dF/dfia is evaluated at /Za=0. 
Comparing with (4.3), the current is Jr = + O(a). Thus we can identify the charge 

density at this order as n = . 

At first order in derivatives and up to 0(a), we have the following invariants: 


{d,, a', 


where d* = and a( = These terms can appear in the effective action, multiplied 

by coefficients which are functions of a*. The number of such terms can be reduced by 
integrating by parts and redefining the relevant coefficients, and the action can be written as 


c(l) _ 


ciT -|- C2/i + cs/ia -|- diT -|- d2H + d^^^EA -|- 0{a^) , 


(4.6) 


where Cj and di are functions of a^, and F^ = is the r-type (physical) electric field. This 

again has the form of (4.3). Combining the contributions from (4.5) and (4.6), the current to 
order 0{aP) is 

+ f , (4.7) 
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with u^jf^ = 0, and 


(4.8a) 

(4.8b) 


= Ffia + {ci,ij,a-C3,T-di)t + (c2,^,-C3,^-(i2)A “ c^d-u, 

f = diAf^'^d^T + d2A^^a^/i + dsAf^'^u^ + diE^^. 

It is reassuring to see the dissipative contributions to the current emerge from the effective 
action formulation. The constitutive relations (4.8) are written in a general “frame”, and 
without using the ideal hydrodynamics equations of motion. The identification of the transport 
coefficients in terms of the parameters of the effective action needs to be done after the same 
expansion in the a-type fields is implemented for the energy-momentum tensor. 

4.2 Neutral fluid 

We now proceed along similar lines in considering the following realization of diffeomorphism 
invariance, x 2?2 in the low energy hydrodynamic regime. Linearized diffeomorphism 

invariance then ensures that the low energy effective action S'eff depends only on the combi¬ 
nation = 2 'Df^ipf, = where is the vector ‘Goldstone-like’ 

mode which transforms by a shift under a-type diffeos, and as a vector under r-type diffeos. 
Again, stands for V”, and the indices can be raised and lowered with the r-type metric. 
The effective action is then 

Ses[gl,i.,glu^V’a] = (4.9) 

It follows that 

Sefi = I \V^Tr{g% - + 0{a^), (4.10) 

where Tr'' = The equation of motion for (/?“ is then V^Tr'^ = 0 as required 

by (3.9a) to 0{a?‘). Again, the unwritten 0{a^) terms must ensure the proper normalization 
of the generating functional in the effective theory. The data available for constructing the 
effective action is given by the set {T, As above, in equilibrium it is more natural 

to determine both the normalized fluid velocity = /3^/y—and the temperature T = 
l/^/^ (where = g^./d^^pn, in terms of a timelike vector /3^. We will not need to make 
assumptions about the metric dependence in the r-sector, so we will work with the conventional 
variables T and u^. 

In terms of this data, there are three scalars at zeroth order in derivatives up to 0{a), 

ai = T , a 4 = ^u, 0:5 = , 

where and = gr'^i^y Thus at zeroth order in derivatives 

= IV^F{T,UQ + 0{a^) = + + (4.11) 

where F is the effective Lagrangian, and the derivatives = dF/d^u and = dF/d^g 
are evaluated at Cu = Cg = 0- We can then read off Tr'^ = 2F^^u^tt^ -|- 2F^^^gr’' -\- 0 (a), or 

rr = (e + p)uV + + 0(a) (4.12) 
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on defining e(T) = and p{T) = 2F^g. Thus at zeroth order in derivatives the 

effective action gives rise to the energy-momentum tensor in ideal hydrodynamics. 

At first order in derivatives and up to 0(a), we have the following invariants: 




where again the dot stands for u^V^. To order 0(a), the coefficients of any structures involving 
derivatives of can only depend on the invariant T, and integration by parts can be used 
to remove all terms with derivatives of from the effective action, so that 


S’!? = jv^r [fiT + f2V^u^ + + 0{a^). 

(4.13) 

The brackets denote symmetrization, -|- etc. To linear order in the 

a-fields, the coefficients /i, /2 are functions of T, and while fs, f 4 , /s are functions of T 
only. Expanding fi, to first order in , we can read off the physical energy-momentum 
tensor from (4.10). Combining the contributions from (4.11) and (4.13), the result to order 
0 {a^) can be written as the standard hydrodynamic decomposition 


rr = £u>^u’' + + (g^u'^ + g^u^) + , (4.14) 

where tt^g^ = 0, = 0, = 0, and + u^u’^. The coefficients are related to 

the parameters of the effective action by 

T = 2Fg„ - 2Fg^ + 2 (/i,^„-/i,g^-/ 3 )T + 2{f2,^^-h,^^)V-u , (4.15a) 

V = + 2 / 1 , 5 / + 2(/2,€,+/4)V-n , (4.15b) 

g^ = /sA'^'S.T + ih-h)u^ , (4.15c) 

= (4.15d) 

Here d is the number of spatial dimensions, and a^’’' = A^"A'^^(Vq,U /3 -|- V^Ua — 2 

is the shear tensor. The derivatives with respect to and are evaluated at = 0. 


Expressions (4.15) should be viewed as constitutive relations in first-order hydrodynamics, 
obtained from the effective action for the hydrodynamic variables and Goldstone fields. The 
energy-momentum tensor (4.14) is a classical 0{a^) quantity and will receive 0{a) fluctuation 
corrections, as expected from (2.17). The classical constitutive relations (4.15) are written in a 
general “frame”, and without using the ideal hydrodynamics equations of motion. The “frame” 
is inherited from the effective action, similar to the thermodynamic frame of Ref. [15]. There 
are two transport coefficients in first-order classical hydrodynamics, the shear viscosity g and 
the bulk viscosity /. They can be identified from frame-invariant tensor = —r/cr^^ 

and scalar 5 = = -/V • u combinations built from the one-derivative terms 

in (4.15) [16]. The identification of the bulk viscosity requires use of the zeroth-order scalar 
equations UudgT^'^ = 0 to relate the structures T and V • u. Making use of these relations, we 
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as well as 


C = -2(/2,«, + ^/4) + 2vUTh,^^+f2,i^-f2,i,) - 2t;s"r(/l,5.-/1,^,-/3) , 


where = dpjde is the speed of sound squared. 


4.3 Charged fluid 

We can generalize this discussion to the case of a charged fluid, by adding a background 
chemical potential. A priori, we can simply add the gauge field data for the r and a sectors 
to the metric data above. However, the additional vector transforms nontrivially under a- 
diffeos, see (3.8b), and thus we need to modify the tensor data appropriately. We will work with 
the following basis of tensor data that has manifest invariance under a-gauge transformations 
and a-diffeos, 

(4.16) 

where X^j, = Cjl ~ (Note that depends only on and is a- and 

r-gauge invariant.) This data is manifestly invariant in the o-sector, and we can proceed to 
build r-gauge and r-diffeo invariants at the appropriate order in the derivative expansion. At 
0{a), the effective action is 

Seff = + + (4.17) 

where as before Tr'^ = and Jr = 6Ses/6A'^. The equations of motion 

for ifa and ipa are 


= 0, (4.18) 

V^r;, + AlVpJP = , (4.19) 


respectively. The conservation of the r-current ensures r-gauge invariance as required. 
There are five scalars at zeroth order in derivatives up to 0{a), 


ai = T, 02 = g, as = ga, 04 = Cu > «5 = , (4.20) 


where ga = u^X% Cu = before. Note that ga is not manifestly r- 

gauge invariant, and we will need to ensure that current conservation in the r-sector imposes 
invariance in the final equations of motion. Thus at this order, 

*5^? = I V^F{T,Cu,(g,g,ga) + --- (4.21) 

which upon comparing with (4.17) gives T/j, = 2F(r,/r), 5 „u^u'^ -|- 2F{T, g)^^^gr'^ -[- 0(a), as 
well as Jr = F(T, + O(a). 

At first order in derivatives and up to 0(a), we can combine the analyses of the previous 
two subsections to find the following invariants; 


{v^x((, 


XaUg, 


U^Xl 


:XaK,V.U^ 
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together with hj = u^d^ai and a'- = XadiJ.(^i- To order 0{a), the coefficients of any structures 
involving derivatives of and can only depend on the invariants T and /r, and integration 
by parts along with redehnitions of the other coefficients can be used to remove these terms 
from the effective action, so that 

'S'eff = j (/iT + C2/i + C3/ia + diT' + 4.2^' + + dAXa^fi 

+ ©(a^). 

(4.22) 

The notation for the coefficient functions has been chosen to match the earlier discussion as 
much as possible. To linear order in the a-helds, the coefficients /i,/ 2 ,C 2 are functions of o;*, 
for z = 1,..., 5, while f^, f^, /s are functions of T and fi only. Following the earlier discussion, 
we can expand all the monomials to 0{a) and read off the energy momentum tensor and 
charge current as follows 

, (4.23) 

= Mu>^ + f . (4.24) 

The coefficients are related to the parameters of the effective action, and for the energy 
momentum tensor are given by 

T = + 2(c2,^„-C2,^g - /6)/i + f 2 ,^g)V-u , (4.25a) 


T" — 2F^g + “^fx^gT + 2c2,gg/i + 2(/2,gg+^/4)V-u , (4.25b) 

q>^ = + {h-h)u^ , (4.25c) 

, (4.25d) 

while the current is 

^ + {fl,^la-C 3 ,T-dl)t + (c2,;,,-C3,^-(i2)/i + (/ 2 ,/., “ C 3 )V-U , (4.25e) 

f = diA^^^'duT + d2Af^’'dufi + dzA^^^'uy + diE^^. (4.25f) 


These constitutive relations are presented in a specific hydrodynamic frame, and it is useful 
to determine the frame-invariant transport coefficients. For the charged fluid, there is one 
tensor, one vector, and one scalar invariant at first order in the expansion [16]. These are 
usually identihed with the shear viscosity r/, the conductivity a and the bulk viscosity C- In 
terms of the coefficients in the effective action, the tensor invariant remains as for the neutral 
fluid and we find again that 

ri = -U, (4.26) 

except that /4 is now a function of both T and 

The scalar invariant for the charged fluid takes the form S = "PT) _ 
in terms of the one-derivative data in (4.25). This invariant depends on the three tensor 
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structures T, -u, but using the two longitudinal ideal hydrodynamic equations = 

F’^PuvJp = 0 and = 0, we can write S = —CV • n, with the coefficient uniquely identified 

with the bulk viscosity. The result is lengthy, so we will not present it explicitly. 

A new feature of the charged fluid is the existence of a vector invariant. In terms of the 
transverse structures in (4.25), it is given by and depends on the four 

transverse vectors tt“, The transverse ideal hydrodynamic equations 

impose only one constraint among these structures. To isolate the charge conductivity as the 
unique transport coefficient in this sector, we require additional constraints that follow, for 
example, from the equilibrium generating functional [14, 15] (or from positivity of the local 
entropy production). Namely, there is in fact only one linear combination of these structures 
which is consistent with the background equilibrium state, A"l^cl^(/r/T) — However, these 
additional constraints are not apparent in the effective action above, which was derived purely 
on the basis of a specific (r, a) symmetry realization. 

5 Discussion 

The top-down construction of the hydrodynamic effective action in Section 4 reproduces several 
features of the generating function for linearized hydrodynamics in Section 2. However, there 
are also certain missing elements, e.g. as noted for charged fluids in the preceding subsection. 
We therefore conclude by listing several questions left open by the present analysis. 

• (1,2) Basis: We wrote down the hydrodynamic effective action by demanding r- and 
a-diffeo and gauge invariance, although only the r-sector symmetry was linearly realized. 
While this is sufficient to reproduce the expected tensor structures in the classical consti¬ 
tutive relations at 0{a), going to higher orders in the a-expansion requires implementing 
the invariance under 1- and 2-sector symmetry transformations, for example in order to 
reproduce the conservation laws (3.9). Manifest (1, 2) gauge and diffeo invariance is not 
straightforward in the (r, a) basis, while the classical hydrodynamic equations are not 
straightforward to represent in the (1,2) basis. Understanding the translation seems 
important for tackling several of the open questions below. 

• Equilibrium constraints: In classical hydrodynamics, the existence of an equilibrium 
state in the presence of sources [14, 15], or (in some cases equivalently) the positivity of 
local entropy production, leads to powerful constraints on the possible thermodynamic 
response and transport coefficients. It is not immediately clear how to think about the 
entropy current from the point of view of the hydrodynamic effective action we have 
described, and in turn the source-dependence of the r-sector fields in equilibrium is not 
manifest. Clarifying these features would, for example, allow the correct identification 
of the electrical conductivity through a constraint among the coefficients of the tensors 

A^^dpn and U". 

• Fluctuation-dissipation constraints: In addition to the equilibrium constraints above, 
which require the coefficients of certain tensor structures to vanish, it is apparent from 
the bottom up construction in Section 2 that there are nontrivial relations between 
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different orders in the a-expansion of the effective action (3.1). For example, the dissi¬ 
pative transport coefficients enter both Jr and /Cr, as required to satisfy the fluctuation- 
dissipation theorem. It is important to understand the origin of these relations, and to 
determine any connection to the vanishing conditions above for specific coefficients in 
the classical constitutive relations. 

• External fields: The charged fluid analyzed in Section 4.3 raises a question of how a- 
sector diffeomorphism invariance can be manifest in the presence of the explicit violation 
induced by the background electromagnetic field. The vector Xa is not manifestly in¬ 
variant under r-sector gauge transformations, and indeed it seems clear that it should 
enter only in the combination = u^Xa- This leaves open the question of how such 
terms should enter at quadratic order in the a-expansion, as is required to restore the 
fluctuation-dissipation relation. 

• Path integral measure: The hydrodynamic effective actions such as (4.21) and (4.22) are 
meant to be used in the path integral with both r-type and a-type dynamical fields. 
Going beyond 0(a) in the effective action requires understanding what the integration 
measure is for the r-type (physical) variables. The microscopic definition implies that 

= 1) s-iid in the low energy regime the nontrivial measure factor required to 
ensure this was (formally) introduced as the functional Ir- For the linearized hydrody¬ 
namics of Section 2 the issue does not arise, as the action is linear in the r-type fields. 
Knowledge of the correct measure is important to determine off-shell interactions of the 
hydrodynamic degrees of freedom. 

• Galilean hydrodynamics: The focus of this paper has been on relativistic hydrodynam¬ 
ics. It would be interesting to develop this approach for systems with Galilean in¬ 
variance, taking advantage of the recent understanding of Newton-Cartan sources in 
non-relativistic fluids [21, 22, 23]. 

Note Added: As this paper was being finalized, the paper [24] appeared on the arXiv. The 
latter work includes a comprehensive classification of non-dissipative transport to all orders 
in the hydrodynamic expansion, and has some overlap with Section 4 of the present paper 
in its treatment of dissipative terms within a GTP-like formalism. It would be interesting to 
understand the relations between these two approaches in more detail. 
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Appendices 


A Shear and sound response functions 

To see that that the generating functional (2.14) does indeed reproduce the required response 
functions given in (2.12), let us go to Fourier space and integrate out the </?“ fields. This will 
give rise to an effective action which is non-local, but real. With the properly normalized 
measure, we have 


For the case at hand, Mij = —7^(k^(5jj—fcjfcj) — ^skikj, hence 


{M-% = - 


7r?k2 

The generating functional now becomes 


1 kikj 


Z[hr,ha\= Dni exp 


IJ 


ki kj 

J 7s k^ k^ 


F,[7r, hn{M-%F,[7r, h^] + 0{h-) 


where 


Fi[7r, h^] = SijTTj - iwuih^i - -wikih^Q 


We have defined 


kiki 


V I 


k2 


kik^ 


Sij — Ari ( 6ij ) -I- As , {S ^)ij — ( 6ij 


1 


kikj'^ _l_ ^ kikj 


k2 


As k2 


where A^ = {—ioj + A* = {—ioo -|- -|- 7sk^)- Setting the sources to zero, we can 

evaluate the rr function (vrjvrj), which is the symmetrized function (half the anti-commutator). 
With the sources set to zero, we have 


(vTiTT^ 


J ■ ■ ■ . 


= / DiTi exp 


1 f .S^M-^S 
— tt' -vr 


TTiTT. 




2 J 2Tw 

The rr two-point function is therefore {iriiTj) = —2Tw , or more explicitly 


(tTjTTj) 


^ kiki\ hkj 


k2 


|As|2 k2 


(A.l) 


This agrees precisely with the hydrodynamic rr function (2.13) found from the ra func¬ 
tion (2.12a). The rr functions involving e can be obtained from the energy conservation 
constraint, 


hh- 

^rr _ ^rr 

« “ a;2 ^i-^j ’ 


k ■ 

/~irr _ ^rr _ /~irr 

I * UJ ^ 


Again, they agree with the hydrodynamic rr functions (2.13) found from (2.12b), (2.12c). 
The rr functions can of course be obtained by varying the generating functional (2.14) with 
respect to the a-type sources, 


= 2i 


. 6^Z[hr, ha 


hh°‘ hh°- 


20 




















As usual, the sources are set to zero after the variation. In order to find the ra and ar functions, 
we need to vary with respect to one r-source and one a-source. Looking at the generating 
functional (2.14), this will bring down one factor of vTj and one factor of dip^. Specihcally, in 
our conventions we have 


Qra _ ■ ^ Z[hr., hg] 

r ^ . 5‘^Z[hr,ha] 

5hl, 


iw (vTi dtip'j) - pdij , 
iwidtpf-Kj) -p6ij . 


(A.2a) 


(A.2b) 


Integrating out ip^ and e in the generating functional (2.14) and setting the sources to zero, 
we have the following correlation functions: 


( 7 ri<p^...}= / DTTkDipf vr^ 99" 


with Mij and Sij dehned above. This can be schematically represented using the combined 
field \a = (vTj, ipl) as 

(AcArf...) = JdX XcXd-- - , 

where the matrix Kab in the 7r(/?“ space is 

A'=(» ), ;,.-i 

\-iS -2TwMj \ iS^-^ 0 j 

The correlation functions are given by (XaXb) = {K~^)ab, so we have 

{TTiTTj) = -2Tw (^S-^MS^-^y , 

{iTiip'^) = i , {p^j) = i > 


The rr function is precisely what we have just evaluated in (A.l) by integrating out the 
field hrst. For the mixed functions, we have from (A.2) 


^ ra 

^TTiTTj 


as well as 


kikj 

«« - 


ioj — 7,jk2 


+ e 


+ 


kikj 


WUJ 


^2 _ ]^ 2^2 _j_ 


^ ar _ 


Sij 


kikj 

“k^ 




—iuj — 7r)k2 


+ e 


+ 


WUj'^ 


kikj 

k^ [a;2 _ k2^2 _ 


P 


The mixed functions agree with the response functions (2.12a) obtained by varying the on- 
shell hydrodynamic equations of motion, including the contact terms. The aa functions all 
vanish due to (pfipj) = 0. 

For the response functions involving the energy density we have 


. 6‘^Z[hr, ha 


(7™ = 2? 

<ria xur 
OUqq oriQ- 


= iw{e dtp ^), G™ = 2i = -iw{^T^ djp^). 

OUoi OUqq 
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The factor of 2 is due to the coupling of /iqo- These can be evaluated by using e = kiTTi/uj, which 
is imposed in our generating functional (2.14), and the answer agrees precisely with (2.12b). 
Similarly, 

era _ »■ /tg] _ ki ^ 

~ A/)’’ ~ / I ^ ’ 

O/loo OtlQQ 00 

where again the factor of 4 is due to the coupling of /ioo- This agrees precisely with (2.12c), 
including the contact term. 
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